Microwave conductivity due to impurity scattering in a d-wave superconductor 
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The self-consistent t-matrix approximation for impurity scattering in unconventional superconduc- 
tors is used to interpret recent measurements of the temperature and frequency dependence of the 
microwave conductivity of YBCO crystals below 20K. In this theory, the conductivity is expressed 
in terms of a frequency dependent single particle self-energy, determined by the impurity scattering 
phase shift which is small for weak (Born) scattering and approaches 7r/2 for unitary scattering. 
Inverting this process, microwave conductivity data are used to extract an effective single-particle 
self-energy and obtain insight into the nature of the operative scattering processes. It is found that 
the effective self-energy is well approximated by a constant plus a linear term in frequency with a 
small positive slope for thermal quasiparticle energies below 20K. Possible physical origins of this 
form of self-energy are discussed. 



Microwave surface resistance measurements on high 
quality YBCO crystals provided some of the early evi- 
dence for unconventional superconductivity in the high 
T c cuprates. The observation of a broad peak in the 
conductivity versus temperature, well below T c [Q, was 
interpreted in terms of a rapid drop in the scattering 
rate due to the disappearance of the inelastic scattering 
which completely dominates the transport near T c . At 
the lowest temperatures, the conductivity is determined 
by scattering due to static disorder. In the context of 
a generalized two-fluid model, the low frequency con- 
ductivity is proportional to the product of a scattering 
time, r, and a density of thermally excited charge carri- 
ers, n n (T). This conductivity was observed to decrease 
roughly linearly with T . Based on this observation, it 
was conjectured Q that r saturated at some large low 
temperature value, implying that n„(T) varies roughly 
linearly with T . Subsequently it was observed Q that 
the penetration depth (and hence the superfluid density 
n s (T) = n s (0) — n n (T)) varies linearly at low T, so that 
the low T behavior of cr(T) and n n (T) are consistent with 
a constant scattering rate. 

However, the conjectured temperature independent 
low temperature scattering time is difficult to understand 
from a theoretical perspective. [Q When vertex correc- 
tions are neglected, the transport scattering rate, 1/r, 
is essentially equal to twice the one-electron self-energy. 
Within the context of the self-consistent t-matrix ap- 
proximation for quasiparticles with line nodes, the one- 
electron self-energy is not expected to be constant at low 
energies, and hence one would expect the low tem- 
perature transport scattering rate to be temperature de- 
pendent. Quite recently, the UBC Group has reported 
measurements of the temperature dependence of the mi- 
crowave conductivity of a high purity YBCO crystal at 



five frequencies in the range of 1 to 75 GHz. || These 
measurements both confirm and extend their earlier re- 
sults. Not only are the observed temperature depen- 
dences consistent with an almost frequency and tem- 
perature independent scattering rate in the impurity- 
dominated regime below 20K, but the frequency depen- 
dence of the conductivity at fixed temperature is consis- 
tent with this same scattering rate. 

The confirmation of the simple picture of an energy in- 
dependent scattering rate is problematic and raises ques- 
tions about our understanding of the nature of excita- 
tions from the high T c superconducting ground state. 
This is particularly so because the crystals under study 
are so clean that one might expect theories based on per- 
turbation theory to apply. Nevertheless we will show 
that the weak energy dependence of the scattering rate 
obtained from the microwave conductivity is inconsis- 
tent with standard simple pictures of scattering of d-wave 
quasiparticles from point impurities, and we will discuss 
various possible physical interpretations of this energy 
dependence. We also suggest that the introduction of 
small quantities of Zn impurities could be used to test 
the predictions of standard "dirty d-wave" theory. 

We begin by comparing the microwave conductivity 
data to the standard model written in terms of the energy 
dependent single particle lifetime. Following the work of 
Hirschfcld et al. H the conductivity may be written as: 
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where Q is the microwave frequency. The pref actor, 
ne 2 /to* A, is related to the temperature dependence of 
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the penetration depth, l/(p X xx (T)) — n s (T)e 2 /m* . 
The integral over Q of cr xx (£l, T) — a xx (fl, 0) is equal 
to l/(no\ 2 xx (0)) - l/(/i A^(T)). Thus the prefactor 
can be determined from the slope of the inverse pen- 
etration depth at low T. The result is ne 2 /m*A = 
-{\rx2/2n )d{l/\l x {T))/dT » W^^wT 1 forthea-axis 
conductivity data of Ref . [|| . 

In order to address the question of whether the data 
can be described within the self-consistent t-matrix ap- 
proximation we need an expression for the function 
F(u>, fi) for a general form of the quasiparticle self en- 
ergy. Given this expression, one can, in principle, invert 
the data to extract the energy-dependence of the scatter- 
ing rate, using Eq. ([]]). We have derived an expression 
for the function F(u>, fi) which is similar to that of Ref. 
[^| . Working in the "node approximation" , expanding the 
quasiparticle dispersion relations around the four d-wave 
nodes, and neglecting the real part of the quasiparticle 
self energy, we find 
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where T(co) is the imaginary part of the quasiparticle self- 
energy. The — ^ limit of this expression is 

F(w,0) = — jr-r tan- 1 -— + -. (3) 

Inserting Eq. (|^) into Eq. (|l|) in the limit of T — ► gives 
/■Km* A, the universal limit [Q, provided that 
the first term in Eq. (||) vanishes when ui — > 0. Assuming 
that F(w) and f2 are both much less than u> gives Eq. (Q) 
below with 1/r = 2r(w). 

fi) is a complicated function of the two fre- 
quencies and of the complex quasiparticle self-energies. 
Hirschfeld et al. B derived the simplified form, 



F{oj, Q) = Im[|w|/(fi - 
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which is a good approximation when the microwave fre- 
quency and the scattering rate are both small compared 
to T and all three energies are small compared to A. 

Surface resistance measurements are converted to con- 
ductivities, using the expression 

J? s (w,T) = ^ w 2 A 3 (a;,T)a 1 (a;,T) (5) 

The temperature dependence of the penetration depth 
at 1 GHz is taken from measurements made at that fre- 
quency. For higher frequencies there is a small frequency- 
dependent contribution to the penetration depth which 



we obtain self-consistently from our fits to the point- 
scattering model. The surface resistance data of Hosseini 
et al. §1 for T < 20 K are reproduced in Fig. |. The 
solid lines in the figure are quadratic fits which allow 
interpolation and extrapolation to T = 0. The T = 
extrapolations of these data, converted to conductivities, 
are shown in Fig. |[ |sj All of the extrapolated data are 
substantially larger than the expected "universal limit" 
M , and it has been suggested H that this discrepancy is 
a real effect due to vertex corrections which enhance the 
universal conductivity at T = 0. In the Born limit, it is 
enhanced by a factor of (t 4i ./t) 2 where r tr is the trans- 
port lifetime. The physical origin of these corrections is 
the fact that, at low T, scattering within a node does 
not change the electric current whereas internode scat- 
tering does. The vertex corrections take the momentum 
dependence of the scattering potential into account by 
correctly weighting scattering among the nodes. 

Experimentally it appears that, in addition to a large 
correction at low T, there is also a frequency dependence 
to the T = limit, with the conductivity exhibiting a 
low frequency peak and a slow roll-off at higher frequen- 
cies. Theoretically, it is unclear what to expect for the 
frequency dependence of the conductivity at T = 0. At 
present there are no calculations available for the fre- 
quency dependence of the vertex corrections. If vertex 
corrections are neglected, then one expects the T = 
conductivity to roll off at a frequency which depends on 
the type of scattering. For unitary scattering, the roll- 
off is at a rather high frequency, some fraction of T c . p"l| ] 
For a frequency-independent scattering rate, 1/r, we find 
that the T — conductivity falls to half its zero frequency 
value at about 2.6/r. Thus, this frequency dependence 
seems to be sensitive to the detailed nature of the scat- 
tering mechanism. 
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FIG. 1. Quadratic fits to the surface resistance data of Hos- 
seini et al. for T < 20 K. 




— i 1 1 

o Extrapolated a (T=0) 
• 1 .2 K Data 



_1_ 



_1_ 



_1_ 



20 40 60 

Frequency / GHz 



80 



FIG. 2. The T = extrapolations of the surface impedance 
data converted to conductivities. Shown for comparison are 
the 1.2 K data. The dashed line indicates the predicted uni- 
versal limit. 

The data for the temperature dependence of the mi- 
crowave conductivity are shown in Fig. [| The data are 
approximately linear in T at all 5 frequencies, although 
there is a small downward curvature in the lowest fre- 
quency data and a comparable upward curvature at the 



highest frequency in addition to the clear non-zero inter- 
cepts. The result of Eq. (|l]) for a frequency indepen- 
dent scattering rate, using the simplified form, Eq. (0) 
for F(u),Q,) given above, is 
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For the low frequency data, a xx at 20 K is about 
50ne 2 /m*A which implies that 1/r is about 0.5 K. 

This semi-quantitative fit to a temperature- 
independent scattering rate is unsatisfying in light of our 
theoretical picture of the scattering of d-wave quasipar- 
ticles. If the scattering is weak, as might be expected 
for crystals exhibiting quasiparticle mean free paths of 
microns, then one would expect a scattering rate lin- 
ear in bj (Born scattering). Indeed the factor of \lo\ in 
F(ui, fl) arises from the same source that would lead to 
a linear scattering rate, namely the linear quasiparticle 
density of states. Alternatively, if the scattering is due 
to a very small density of unitary scatterers, then theory 
predicts a scattering rate proportional to l/u>. Of course, 
there is always the somewhat unnatural possibility that 
some intermediate type of scattering might account for 
the data. The self-consistent equation for the quasiparti- 
cle self-energy depends on the density of impurities and 
on their phase shift 5, where S « 1 corresponds to Born 
scattering and 5 w 7r/2 to unitary scattering. We have 
examined the possibility that an intermediate value of the 
phase shift could give rise to an effectively energy inde- 
pendent scattering rate. What we find is that, regardless 
of the value of S, a quasiparticle self-energy with a mag- 
nitude on the order of 10~ 3 A will have an energy depen- 
dence over the relevant range of energies, < u) < A/10, 
which is inconsistent with a constant scattering rate. 

Since the data show noticeable curvature and hence 
do not follow Eq. (||) perfectly, we have fit the tempera- 
ture dependent part of the microwave conductivity using 
a more general, frequency dependent form of T(ui). We 
model T(u>) as a linear combination of weak (Born) and 
strong (unitary) scattering, and we also allow for the pos- 
sibility of a frequency independent (Drude) component. 
Thus r(w) is written as the sum of a linear term, r^w/A, 
plus a term of the form T u A/\/u) 2 + r„A which mimics 
the effect of unitary scattering, plus a constant, To. We 
found that it made little difference in the fits whether 
Eq. (||) or (^) was used, and so we worked with the sim- 
pler expression, Eq. (Q), to fit the temperature-dependent 
part of the conductivity (with cr(fi, T = 0) subtracted 
off). The best global fits, in which data at all frequencies 
and temperatures were fit to a single set of parameters, 
are shown by the solid lines in Fig. ([|). The fits provide 
a reasonable model of the evolution of the curvature of 
a vs. T as a function microwave frequency. The values 
of the parameters which gave the best global fit were: 
T u = 0, T D = 0.33 K, and T B = 0.69 K. 
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FIG. 3. The best global fit to the temperature-dependent 
part of the conductivity below 20 K. 

In addition we performed individual fits to the temper- 
ature dependence of the conductivity at each frequency 
with our model. The results for the best-fit parameters 
are shown in Table (||). Only the 75GHz data were com- 
patible with a non-zero value of the unitary scattering 
parameter, T u and this fit also required a large value 
of the linear parameter, T b ■ The frequency dependence 
of the one-particle self-energy, corresponding to each of 
these fits, is shown in Fig. (||). We conclude that the 
one-particle self-energy consists of a constant term plus 
a small but clearly positive slope. The most notable fea- 
tures of this conclusion are the apparent absence of a 
resonant peak centered at zero frequency that one would 
expect from unitary scatterers, and the increase of the 
self-energy with increasing frequency. 
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TABLE I. Best fit parameters to individual data sets using 
the point scattering model. 



FIG. 4. The frequency dependence of the one-particle 
self-energy, as obtained from fits to individual data sets. The 
data is best described by a constant + linear self-energy. 

It is relevant to review the common rationale for ex- 
pecting unitary scattering to dominate the transport 
properties in YBCO. jl2| Early measurements on both 
crystals and films generally exhibited a relatively flat 
temperature dependence at low T, which reinforced the 
then commonly held belief that the gap had an s-wave 
symmetry, although the temperature dependence was 
more quadratic than exponential. (l^J^] The observation 
of a linear temperature dependence of the penetration 
depth by Hardy and co-workers || in very clean samples 
of optimally doped YBCO, which provided strong evi- 
dence for a d-wave gap, also demanded some consistent 
explanation of the earlier results. Hirschfeld and Golden- 
feld Jljj pointed out that, if the extra impurities in ear- 
lier samples were unitary scatterers, then a d-wave gap 
was consistent with the approximately quadratic temper- 
ature dependence observed at low T in these materials. 
By comparison, the amount of Born scattering required 
to give the observed range of quadratic temperature de- 
pendence would also be expected to lower T c drastically. 
They also argued that unitary scattering accounted for 
the observed linear dependence of the low T thermal con- 
ductivity in dirty samples. 

Although the arguments given above account for the 
behavior of dirty samples, i.e. samples with a significant 
temperature range over which the penetration depth is 
quadratic in T, it is not at all clear that they should ap- 
ply to the samples used in these experiments, which are 
extremely pure and which exhibit a linear temperature 
dependence of the penetration depth down to the low- 
est temperature studied. As we have seen, the transport 
scattering rate in these materials is less than 0.5 K at low 
T, and the contribution obtained from extrapolating the 
linear term in the scattering rate, from the fits above, to 
T c would be less than an additional Q.5K. These rates are 
completely negligible compared to the inelastic scatter- 
ing rate at T c which is roughly 100 K. It seems inappro- 
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priate to attribute this small amount of scattering to a 
single mechanism in these samples where the low T quasi- 
particle mean free path is several microns. More likely 
the scattering is due to a combination of mechanisms, 
including rare strong-scattering by point defects, some 
more ubiquitous weak long-range Born scattering, due 
to oxygen disorder and weakly scattering lattice defects, 
and scattering by extended defects such as remnant twin 
boundaries. Even different kinds of point impurities may 
have distinctly different kinds of scattering properties. 
This is clear from studies of deliberate impurity doping, 
where a single type of impurity dominates the transport. 
Work in this area has shown that although Zn appears 
to behave as a unitary scatterer, other impurities such 
as Ni and Ca seem to be much weaker scatterers. 
Whatever the mechanism or combination of mechanisms 
in the high purity samples, it appears to result in an al- 
most frequency independent quasiparticle lifetime in the 
temperature range which has thusfar been probed. 

Hcttler and Hirschfeld (HH) have recently proposed 
ijlOlfLif that the apparently frequency independent part 
of the self-energy is the result of a resonance that arises 
away from u — due to suppression of the energy gap 
around the impurity. This finite frequency resonance is 
superimposed on the zero frequency resonance normally 
associated with unitary scatterers in d-wave supercon- 
ductors. They associate the gap relaxation resonance 
with the frequency-independent self-energy that we infer 
from the data. We agree that gap relaxation must occur 
around impurities, and that it will affect the quasiparti- 
cle self-energy. It is less clear that it will give rise to a 
second resonance. jl6| In any case the effect of the zero 
frequency resonance which remains in HH's calculation 
would be to suppress the conductivity at low T in a way 
that is somewhat inconsistent with our data. Specifically, 
the lowest temperature (T < 5 K) data at 1 GHz have 
the opposite curvature to that of HITs fit. 

In conclusion, we find that, in the simple picture in 
which the microwave conductivity is determined by the 
imaginary part of the quasiparticle self-energy, the self- 
energy for low-energy quasiparticles acts roughly like a 
constant plus a small linear term with positive slope, and 
this behavior is inconsistent with existing theories of a 
single kind of point-impurity scattering. Below 20 K, the 
constant is less than 0.5 K. We attribute this behavior to 
the cumulative effect of scattering by a variety of dilute 
and/or weak scatterers. The way to test this hypothesis 
is to introduce a single type of controlled disorder, such 
as Ni or Zn impurities. These added impurities will cause 
R s to decrease and to exhibit a frequency and temper- 
ature dependence characteristic of the added impurity. 
The fact that the starting materials are extremely pure, 
means that the doping required to do this can be very 
small, so that the experiments will probe the effects of 
well-separated dopant atoms. Such experiments may also 
provide further insight into the nature of vertex correc- 



tions which depend on the momentum dependence of the 
impurity potential. 
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